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In recent years ihe semi-Dirac semi-Weyl semi-metal has been of interest due to its naturally 
occurring point Fermi surface and the exotic anisotropic band-structure near the Fermi surface, 
which is linear (graphene-like) in one direction of the Brillouin zone, but quadratic in a direction 
perpendicular to it. In this paper the effect of a magnetic adatom impurity in a semi-Dirac system 
is studied. As in a metal, the magnetic impurity in a semi-Dirac system interacts with the sea of 
conduction electrons and gives rise to magnetism. The transition of the semi-Dirac system from the 
non-magnetic to the magnetic phase is studied as a function of the impurity energy, the strength of 
hybridization between the impurity and the bath as well as that of the electron electron interaction 
at the impurity atom. The results are compared and contrasted with those of graphene and ordinary 
metal. Since the semi-Dirac and the Dirac dispersion share similar features,e.g, both are particle hole 
symmetric and linear in one direction, the two systems share resemblances in their characteristics 
in the presence of a magnetic impurity. But some features are unique to the semi-Dirac dispersion. 


I. INTRODUCTION 

It has been about about a decade since a graphene 
layer has been isolated in the laboratory successfully. 
Its unique linear (massless Dirac, properly called 
Weyl)low ener gy band structure has been studied 
extensively.13^ One of the low energy features of 
graphene is the point Fermi surface aspect. In the 
graphene bandstructure, there is gap thoroughout 
the Brillouin zone except for at the point where the 
bands touch each other. The Fermi (E^) lies ex¬ 
actly at that point. Away from the Fermi surface 
the bands open up in a linear manner for as long 
as the energy is not too large. Another new point 
Fermi surface system, named ‘semi-Dirac’ or ‘semi- 
Weyl’ has been of interest in recent years. It has 
the aspects of both the conventional and the un¬ 
conventional dispersions. The dispersion is linear 
(“massless”, Dirac-Weyl) in one of the directions 
of the two-dimensional (2D) layer, and is conven¬ 
tional quadratic (“massive”) in the perpendicular 
direction. At directions between the axes the disper¬ 
sion is intermediate and highly direction-dependent. 
Volovik obtained such a spectrum at the point node 
in the A-phase of superfluid ^He [3] and studied its 
topological robustness.13 Montambaux’s group dis¬ 
covered the semi-Dirac spectrum in a model based 
upon the honeycomb lattice geometry, which also 
happens to be crystal geometry for graphene.^ The 
model has a broken rotational symmetry such that 
hopping to two nearest neighbors is t but to the third 
neighbor is t'. When f' differs from t, the graphene 
“Dirac points” move away from the K and K' points, 
and for the specific choice t' = 2t they merge, result¬ 
ing in the semi-Dirac spectrum. This group began 
a study of low energy properties of such a system,!^ 


which was continued by Banerjee et alW. The point 
Fermi surface along with the anisotropic dispersion 
(massive in one direction, massless in the orthog¬ 
onal direction, and hybrid of the two in between) 
characterizes a semi-Dirac dispersion in an unique 
manner. The semi-Dirac dispersion has recently 
been experimentally realized in the context of optical 
latticed. There ultracold atoms, which mimic the 
electrons in a real condensed matter sustem, have 
been studied in an adjustable honeycomb lattice en¬ 
vironment. By tuning the hopping parameters to the 
right anisotropy the merging of the Dirac points into 
a single semi-Dirac point has been achieved. From 
the simulation standpoint Pardo and PickettP^ re¬ 
ported the finding of the semi-Dirac dispersion in 
ultrathin (001) VO 2 layers embedded in Ti 02 using 
Density functional theory calculations. 

In this paper we investigate how a single magnetic 
impurity hybridizes with the conduction electrons 
in a semi-Dirac system to give rise to the magnetic 
non-magnetic phase transition. This study has pre¬ 
viously been conducted in the context of an ordinary 
metaP^ and graphen^i^l. A metal with a magnetic 
impurity in it does not show any magnetism if the 
energy of the impurity atom is larger than the Fermi 
energy. But a somewhat unusual non-magnetic to 
magnetic phase transition takes place when a mag¬ 
netic impurity is embedded in the graphene crys¬ 
tal. Due to the anomalous broadening of the impu¬ 
rity density of states, magnetism can be observed in 
graphene even when the energy of the impurity atom 
is larger than the Fermi energy. In this paper it has 
been investigated how this issue as well as various 
other ones play out in the context of the semi-Dirac 
system with a magnetic impurity in it. The com¬ 
plete Hamiltonian for this problem consists of a 2 


band terirP representing the semi-Dirac conduction 
electrons, the onsite energy term of the localized im¬ 
purity state, the term representing the hybridization 
of the impurity with the conduction electrons and fi¬ 
nally the term representing the Coulomb interaction 
between the electrons at the impurity site. In this 
paper it is investigated how the semi-Dirac system 
makes a transition from a magnetic state to a non¬ 
magnetic one and vice versa as parameters like the 
impurity energy eo, the strength of the electron elec¬ 
tron Coulomb interaction U, and the the impurity 
conduction electron hybridization parameter V are 
allowed to scan the parameter space. 

II. THE SEMI-DIRAC SYSTEM WITH 
IMPURITY 

In this section we will discuss the four major com¬ 
ponents, as mentioned in the previous section, of 
the Hamiltonian for the semi-Dirac system with a 
magnetic impurity in details. Those are denoted by 
Htb, Hv, Hf, and Hu respectively. Mathematical 
expressions for each of the above terms are given 
below. Htb, the two band tight binding part of the 
Hamiltonian, which corresponds to the semi-Dirac 
conduction electrons, is written as 

Htb = E Vi 2 {k)c\,^{k)c 2 ,aik) + h.c, (I) 

k,<7 

where c^ 2 ) a^k), and Ci(^ 2 ),a{k) are the creation and 
the destruction operators of the conduction electron 
with the spin index cr; I and 2 being the band in¬ 
dices. The off-diagonal element Vi 2 {k) is given by 

P 

Vi2{.k)=vk,-i^, ( 2 ) 

Zm 

V, and m are the characteristic velocity and the mass 
parameters of a semi-Dirac system {h has been set to 
1). Eq .[^is the second quantized, spin added version 
of the single particle, spinless semi-Dirac Hamilto¬ 
nian appearing in [5]. The part of the Hamiltonian 
that represents the hybridization between the con¬ 
duction electrons and the impurity energy state, de¬ 
noted as Hv, is given bjnn 

Hv = + /],feC 2 .<x(fc) + h.c), (3) 

^ k,(T 

where fl j,, and fa,k are the creation and the destruc¬ 
tion operators corresponding to the impurity state 


of spin a. N denotes the normalization factor The 
onsite energy term for the impurity is given by 

Hf = eoJ2flkU,k. (4) 

k,(7 

Co is the energy of the localized impurity state. So 
far although the spin indices appear in the Hamil¬ 
tonian, operators with opposite spins are decoupled 
from each other. The introduction of the Coulomb 
interaction term Hu, the real space expression of 
which is given in the following changes that. 

Hu = Uflf^flfi, (5) 

It is clearly seen from the above equation that oper¬ 
ators with opposite spins interact with each other. 
U appearing in the expression for Hu is a constant, 
which denotes the strength of the Coulomb repul¬ 
sion between the electrons of opposite spins at the 
impurity site. The complete Hamiltonian (H)is ob¬ 
tained by adding together all the individual pieces 
{Htb, Hv,Hf, and Hu), that have been introduced 
so far. 

H = Htb +Hv + Hf + Hu (6) 

Due to the presence of the Coulomb interaction term 
Hu, the Hamiltonian in Eq. can’t be solved ex¬ 
actly. In order to render Hu more tractable we re¬ 
sort to the mean field technique. The mean field 
decomposition of Hu is given as follows. 

Hu « mflMflU) + {flf^)flU - (/tVT)(/lA)](7) 

The last term in the above equation is a constant 
and can safely be ignored. The first two terms have 
a common trait: both of the them have the number 
operator of the impurity state of a given type of spin, 
multiplied by a scalar which is the average of the 
number operator of the opposite spin. Such a term 
when grouped with the on-site energy term for the 
same spin specie appearing in Hf produces a scalar 
like Co- defined as: = eo + Uu-a, where Ua- denotes 
the average of the number operator of spin a. —a 
indicates the spin type opposite to that indicated 
by cr. Hence under the mean field approximation, 
adding Eq. and Eq. together we obtain 

H'f = Y.^.lljTk- ( 8 ) 

k,(7 

The final mean field expression of the Hamiltonian 
is obtained by adding together Eqs. B 13 and [^as 
follows 

Hmf = Htb + Hv + H'f, (9) 
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Instead of using the exact Hamiltonian given by 
Eq. 1^ the simplified mean field version of it as given 
by the above equation will be used for all the calcu¬ 
lations in this paper. 


III. GREEN’S FUNCTION FOR THE 
IMPURITY STATE 


|e| < D. It turns out that the integral can be evalu¬ 
ated in closed form. From Eqs.[TT]and [T^one finally 
obtains 
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,( 13 ) 


where /(w) is given by 


In order to find out if the semi-Dirac system is 
in the magnetic phase or not for a certain parame¬ 
ter set, one needs to solve for impurity occupation 
numbers corresponding to two different spins in a 
self consistent manner. If the difference between the 
occupation numbers is non-zero that would imply 
that the system is magnetic phase, whereas a zero 
difference between them would be indicative of a 
non-magnetic phase. In order to determine impu¬ 
rity occupation numbers finding the Green’s func¬ 
tion for the impurity state is essential. From the 
Green’s function the impurity density of states can 
be obtained, which upon integration will give an ex¬ 
pression for the impurity occupation number. The 
impurity Green’s function ^ for the semi-Dirac 
system is defined as Gj^ g. =< T[faf^] >,T denoting 
the time ordering operator. An expression for G^j ^ 
derived using the above definition is given bjE^ 


_ 


w - e^r - ^ - iOt ■ 


( 10 ) 


^ in the above equation is the self energy is writ¬ 
ten in terms of non-interacting Green’s function as 
follows 
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( 11 ) 


The non-interacting Green’s function G^ g. in the 
above equation is given by 
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At this point a brief discussion about the choice 
of D is warranted. D is ascertained following the 
Debye prescriptiorPl^. According to the Debye pre¬ 
scription, D is determined based on the conservation 
of number of states in the Brillouin zone. D for the 
semi-Dirac system is a function of the adjustable pa¬ 
rameters m and v. Since m and v need to be deter¬ 
mined empirically from the real physical system, D 
can differ from one semi-Dirac system to another. 
To be able to compare the results with graphene, 
the value of D for the semi-Dirac is chosen to be the 
same as that in graphene, for which D ^ 1 eV^. 
[The semi-Dirac parameters m and v depend on the 
tightbinding parameter t, and the bond-length cP. 
Using the same numerical values for t and a as are 
used for graphene, D for semi-Dirac has been veri¬ 
fied to be of the same order of magnitude as that in 
case of graphene. Hence the aforementioned choice 
of D is reasonable.] As long as the Fermi energy /r is 
much less compared to D, the physics does not get 
affected by a specific choice of the value of D. 

The localized density of states for the impurity 
atom is obtained from the imaginary part of the 
Green’s functiorPl^. 


From Eqs. fldjand [^we obtain 

3n y’’ 


(15) 
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PffA^) = - 




7r(a.Z-i(a;)-e.)2 + (^^)2 
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where is the square of the quasi-particle energy 
of the semi-Dirac system, which in the momentum 

space takes the form £^ = A + The summa¬ 
tion over momentum in Eq. El is converted into an 
integral involving the density of states of the semi- 
Dirac quasi-particle. The integration is carried out 
from —D to D, where D is the cut off energy. The 
Semi-Dirac dispersion is assumed to remain valid for 


where Z ^(w) is given by 

qT/2 

= i + (17) 

The impurity density of states appearing in 
Eq. looks very different when compared to the 
Lorentzian impurity density of states for an ordinary 
metaP^. It also differs from the impurity density of 
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states of graphene, e.g, in the presence of a square 
root of the absolute value of the energy uj instead 
of simply the absolute value of oj in the numerator. 
The occupation {n^) of the impurity is given as the 
integral of p//,dr up to the Fermi energy as follows 

Tier = f dujpff{uj), (18) 

J —OO 


The above equation is written as the sum of two sep¬ 
arate integrals over the energy ranges (—00 to 0) and 
(0 to p) respectively. Under the assumption that the 
semi-Dirac dispersion is valid from —D to D, as was 
explained before, in one of the integrals —00 should 
be replaced by —D. p being much smaller than D, 
the upper limit of the integral remains unchanged, 
from Eqs. [T6|and [TSl we obtain 

ria = Ii -\- 12 , (19) 


where 
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where A = (37r/4)(U^/Zl5). rio- appearing in the 
left side of Eq. [I^ depends on n_cr through the term 
Ccr appearing on the right side of the same equation. 
Hence Eq.p^needs to be solved self consistently. For 
a given uj the integrands in the above equati ons are 
large when Z~^(lu)uj « Co-. Referring to Eq. 17 and 


keeping only the dominant term in the expression of 
Z~^{uj)uj for a small A and not too small uj, it is 
clear that uj ^ Ccr- Hence Z~^{uj) can be approxi¬ 
mated by Z~^ = Z~^{ea-). With this approximation 
h and I 2 can be evaluated in closed forms. It has 
been checked that the results obtained with the ap¬ 
proximation match rather well with those using the 
complete uj dependence of Z~^(uj). To keep the fi¬ 
nal closed form expressions of the integrals neat, the 
following quantities are defined. 


a = A -|- 2Z^ €a, 


(21a) 


6 = A" - 2Z-ie,. 


(21b) 


Evaluating Ii and I 2 as given by Eq. 20 we obtain 
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L[ = 2Z-2(g^) [^+(23g) 


L '2 = (23h) 

Her is obtained self consistently from Eq. [fusing the 
Matlab minimization routine ‘fsolve’. A shift in the 
difference between the values of na for two different 
cr’s corresponding to two different spins, from the 
zero to a non-zero value implies that the the system 
has made a transition from the non-magnetic to the 
magnetic phase. Defining the variables x, and y as 
X = and y = (p — to)/U, the boundary (which 
henceforth will simply be referred to as the transi¬ 
tion curve) separating the magnetic phase from the 
non-magnetic one is plotted as a function of x and 
y [details to be given in the ‘Results’ section]. The 
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more the number of grid-points in the x — y plane 
used to produce the plots, the sharper is the transi¬ 
tion curve, but that happens at the cost of somewhat 
increased program runtime. Hence an optimal num¬ 
ber of grid points have been used to get a fairly good 
idea of how the transition curve should look like for 
a particular parameter-set. Then smoothing algo¬ 
rithm has been employed to smooth out any small 
irregularities that might be present at the boundary 
region. But no smoothing techniques were necessary 
to produce the plots for the occupation numbers for 
different spins as functions of /i. Both positive as 
well as negative values of eo have been considered 
for the plots. It turns out for the negative values 
of eo the convergence of the self consistent solution 
using fsolve optimization routine can be somewhat 
poor if one works with the closed form approxima¬ 
tions (as given by Eq. 22) of the integrals appear¬ 
ing in Eq. A possible reason for that is dur¬ 
ing the optimization process Sa in the expression for 
Z~^ can go to zero resulting in making Z~^ a large 
number and thereby affecting the convergence in a 
non-physical way. This does not happen when £o is 
positive, since all the terms in the expression of 
are positive and the possibility of £a going to zero 
does not arise. Hence for the negative values of Eq, 
an un-approximated direct numerical integration of 
Eq.ll^has been resorted to, for producing the plots. 


IV. RESULTS 

In Figj^the transition curves are plotted as func¬ 
tions of X and y for positive as well as negative val¬ 
ues of cq. As in the Dirac dispersion, they are not 
symmetric about y = .5. There is also asymmetry 
between eo being positive and negative. The fact 
that both the semi-Dirac dispersion and the Dirac 
dispersion are not symmetric about y = .5 can be 
attributed to the particle-hole symmetry breaking 
due to presence of the impurity. From the plot in 
the top of Figj^t is observed for the positive value 
of eo, the magnetic region barely crosses the line 
y = 0, indicating that the system is magnetic even 
when the impurity energy is slightly larger than the 
Fermi energy. For an ordinary metal, the transition 
curve never goes below y = 0. In case of Dirac, 
the magnetic region extends much further into the 
region where y is negative with an overall negative 
slope, which is in direct contrast to the semi-Dirac 
case. When y is negative, the on-site impurity en¬ 
ergy is larger than the Fermi energy, but because of 
the anomalous broadening of the impurity density of 



Semi-Dirac magnetic/non-magnetic boundary plot for the parameter set: 

D = 7 0, X„^er0.0, x,ppe,=10.0, y|,^g,=-0.3, y =1 -f, X step size=0.01; y step siz 
(V=1.0, epsQ=-0.2), (V=1.0, epsQ=-0.3), (V=0.2, epsQ=-0.2), 
red corresponds to the first set, blue and green to the second and the 3rd sets respc 



FIG. 1: Magnetic to non-magnetic transition curves for 
the following parameter sets: Top figure: red corre¬ 
sponds to £o = .029D, V = .14D; blue corresponds 
to £o = .0431), V = .14D, and green corresponds to 
£o = .029D, V = .03D. Bottom figure: red corre¬ 
sponds to £o = —.029D,V = .14D; blue corresponds 
to £o = —.0431), V = .141), and green corresponds to 
£0 = -.029D,V = .03D 


states magnetism is still manifested. For the semi- 
Dirac the confinement of the magnetic region pri¬ 
marily to the positive values of y can be attributed 
to the anomalous but somewhat weak broadening 
and the resulting lack of ease with which the impu¬ 
rity density of states crosses the Fermi energy, com¬ 
pared to that of the Dirac case as is explained in the 
following. The impurity density of states of neither 
the semi-Dirac nor the Dirac dispersion is a simple 
Lorentzian. For the semi-Dirac it falls off as , 

as can be seen from Eq. whereas for the Dirac 
system, the tail of the impurity density of states goes 
as So, for the semi-Dirac the impurity density of 
states decreases faster, and hence crosses the Fermi 
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energy with more difficulty, compared to the Dirac 
case. But the semi-Dirac impurity density of states 
does not fall off as fast as a metal, for which the 
density of states being Lorentzian, falls off as 
That explains the difference in the behavior of the 
transition curves in the region of negative y close to 
0 between the semi-Dirac and the metal. Whereas 
y = 0 corresponds to eg being equal to fx, y = 1 
corresponds to the case when the impurity state is 
doubly occupied and the larger of the two impurity 
energy levels is equal to the Fermi energy. From the 
plot in the bottom of FigJ^ it is observed that the 
transition curve for the semi-Dirac dispersion for the 
negative impurity energy, although somewhat simi¬ 
lar in appearance to the Dirac, barely crosses the 
y = 1 line(it extends considerably beyond y = 1 for 
the Dirac case). For small x, the transition curve 
shows a kink, which implies that the curve near 
cc = 0 goes downward like the ordinary metal or 
the Dirac. This universal behavior in the transi¬ 
tion curve is brought about by the irrelevance of the 
specific nature of energy momentum dispersion of 
the conduction electrons in the small x limit, since 
a small value of x indicates weak coupling of the 
magnetic impurity with the conduction electron. A 
small value of x corresponds to large U. Hence it 
is concluded that a large U causes the system to 
behave universally. It turns out that compared to 
an ordinary metal, a much smaller value of U can 
send the semi-Dirac system to magnetic state. For 
V = \, and eg = .3 the magnetism sets in when U 
is as low as .35 (This is determined as follows. For 
a given V, and D, delta is calculated. Then from 
the top plot in Figj^the maximum value of x is as¬ 
certained for which magnetism still exists, x being 
a function of U, and everything else in the expres¬ 
sion of X like D and A being known, the minimum 
U is thus determined from the maximum possible 
value of X.) The minimum U thus obtained for the 
semi-Dirac turns out to be somewhat larger than the 
critical U of graphene for the same parameter selP^, 
but much less compared to U in an ordinary tran¬ 
sition metal. Hence transition material atoms, not 
magnetic in isolation, can become magnetic when 
introduced in a semi-Dirac system. Next the oc¬ 
cupations of the impurity energy level for different 
types of spins are plotted as functions of Fermi en¬ 
ergy. In Figthe impurity occupations are plotted 
for for two different values of a;’s: 2.1 and .5. Both 
the top and the bottom plots correspond to the the 
following parameter set: eo = .2, V = l,andD = 7. 
X = 2.1 for the top plot and x = .5 for the bot¬ 
tom one. For x = 2.1 it is observed that there is a 


Semi-Dirac n . versus mu plot for the parameter set: 
D=7.!)° V^1.0, epSj=0.2, x=2.1 U=0.4 



Semi-Dirac versus mu plot for the parameter set: 
D=7.0^ V=1.0, epsQ=0.2, x=0.5 U=1.8 



FIG. 2: tier versus y plot. The top and the bottom figures 
have identical parameter set except for the values of x. 


magnetic to non-magnetic transition above y = .2. 
This is different when compared to the Dirac dis¬ 
persion, for in the Dirac dispersion the transition 
happens at y somewhat lesser than .2 for the same 
parameter selfS^. There are also other differences, 
like the range of y over which there is a non-zero 
magnetic moment is twice for the semi-Dirac system 
than that in Dirac. Also, for semi-Dirac the mag¬ 
netic moment is slightly larger (~ .GS/is) compared 
to Dirac’s ~ .byj^or the same parameter set. The 
larger range of magnetism for the semi-Dirac can be 
attributed to its larger value of U due to the ap¬ 
pearance of (instead of as is the case for 
Dirac)in the denominator of the expression for A. 
That means for a given set of V and D, A for the 
semi-Dirac is larger. Hence for a given value of x, 
U for the semi-Dirac is also larger as can be seen 
from the expression x = A large U is more 

likely to favor magnetism, which can be seen con¬ 
sidering the extreme case of t/ —> oo, when due to 
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the large energy cost of double occupation, the im¬ 
purity site will always be occupied by a single spin 
resulting in a non-zero magnetic moment. In the 
same figure, at the bottom the occupation numbers 
are plotted for a: = .5. cq, as well as all other pa¬ 
rameters remain the same as in the top. Compared 
to the figure in the top, the non-zero magnetic mo¬ 
ment is over a larger range of /r and the magnetic to 
non-magnetic transition is rather abrupt for the bot¬ 
tom figure. The magnetic moment in the magnetic 
region gets as large as All these features are 

rather similar to the Dirac case, except for a some¬ 
what larger range of ^ for which the system stays in 
the magnetic phase. When eg is negative, the tran¬ 
sition takes place for a value of ^ which is barely 
negative (not shown in the figure). In that respect 
the semi-Dirac system is different from Dirac; for in 
Dirac, the non-magnetic to magnetic transition takes 
place at a larger negative value of fi. The magnetic 
susceptibility for the impurity is computed using the 
following expressiorf^ 


x = -AY. 


G 


dUa- 


1 - U 


dri—a 

de_„ 


df-a 1 — drbi 


(24) 


where is the Bohr magneton. It been checked 
that the susceptibilities (not shown in the plot) show 
sharp peaks for those values of fj,, at which the non¬ 
magnetic to magnetic transition takes place. 


tern. The magnetic impurity couples with the bath 
of electrons of the anisotropic semi-Dirac system to 
give rise to a unique magnetic to non-magnetic tran¬ 
sition phenomenon. The unconventional energy de¬ 
pendence of the density of states of the semi-Dirac 
conduction electrons results in a rather unique hy¬ 
bridization with the magnetic impurity atom, mak¬ 
ing it very different from either graphene or an or¬ 
dinary metal. In this paper a specific model Hamil¬ 
tonian for the semi-Dirac system has been used. A 
future direction will be to use other types of Hamil¬ 
tonian for the semi-Dirac system and check the effect 
of magnetic impurity for that system. Semi-Dirac 
is a relatively new discovery, which exists in theory 
and in simulation. In coming years as understand¬ 
ing about this unique system gets better, it will be 
clearer if the Fermi energy can be controlled suc¬ 
cessfully by the application of some kind of a gate 
voltage, a direct application of which will render the 
semi-Dirac material valuable for spintronics. This 
paper is a first attempt to understand the interac¬ 
tion between a magnetic impurity and a specific type 
of semi-Dirac system, hoping that it will give an idea 
about the magnetic to nonmagnetic phase boundary 
and magnetic moment formation dependent on the 
parameters like impurity energy, impurity to con¬ 
duction electron coupling, and Fermi energy of the 
system. 


V. CONCLUSION 

In this paper we have discussed the Anderson’s 
impurity problem in the context of a semi-Dirac sys- 
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